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Abstract. Non-abelian energy loss in quark gluon plasmas is shown to lead to 
novel hadron ratio suppression patterns in ultrarelativistic nuclear collisions. 
Here we investigate pion and kaon production in pp and AA collisions in a 
perturbative QCD frame, suppression pattern and hadron ratios. The K~ / ' K + 
and K + /ir + ratios are found to be most sensitive to the opacity (density) of the 
plasma. Experimental data indicate that the fragmentation dominated pQCD 
region will be reached only at higher px', in an intermediate px region other 
particle production mechanisms dominate the K/t: ratios. 
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Energy loss of high energy quark and gluon jets penetrating dense matter pro- 
duced in ultrarelativistic heavy ion collisions leads to jet quenching and thus probes 
the quark-gluon plasma formed in those reactions [ [H 13 [3l [4], [5l [6J, [TJ [8]. The 
non-abelian radiative energy loss suppresses the hadron production yield in the mo- 
mentum range 2-3 GeV/c < p T < 15-20 GeV/c at RHIC energies HP]. 
Since quark and gluon jets suffer different energy losses proportional to their color 
Casimir factors (4/3, 3), the energy loss effect could be investigated through the jet 
fragmentation pattern into hadrons with different flavor content. Here we present 
a perturbative QCD based calculation of different hadron ratios in p + p collisions [ 
111! 112] . then extend our calculation for heavy ion collisions and include jet energy 
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loss. 

Let us summarize our basic knowledge about non-abelian energy loss in hot 
dense matter. First estimates [[TJ [5] suggested a linear dependence on the plasma 
thickness, L, namely AE 1 — 2 GeV(L/fm), as in abclian electrodynamics. In 
BDMS [ [3 |4] , however, non-abelian (radiated gluon final state interaction) effects 
were shown to lead to a quadratic dependence on L with a larger magnitude of AE. 
Similar results have been obtained from light-cone path integral formalism [0 [6|. 

In the GLV formalism, applying opacity series (see Refs. [[TJ [5] ) , finite kinematic 
constraints were found to reduce greatly the energy loss at moderate jet energies. 
On the other hand, the quadratic dependence on L has been recovered in wide en- 
ergy region. As we show below, the kinematically suppressed GLV energy loss turns 
out to depend approximately linearly on the jet energy, E. This linear dependence, 
AE oc E, unfortunately leads to a very weak p^-dependent suppression of the trans- 
verse momentum distributions for the kinematic range accessible experimentally at 
RHIC. Our focus is to investigate whether the pr dependence of the particle ratios 
may help to pin down the jet quenching mechanisms. Our primary candidates are 
the measurable K + /ir + and K~ /K + ratios. 

Non-abelian energy loss in pQCD has been calculated analytically in two limits. 
In the "thick plasma" limit, the mean number of jet scatterings, n — L/A, is assumed 
to be much greater than one. For asymptotic jet energies the eikonal approximation 
applies and the resummed energy loss (ignoring kinematic constraints) reduces to 
the following simple form [OH HI EJ [5]: 

a T-i C R a s LV . 

AEbdms = — -; r — v , (1) 

4 A g 

where Cr is the color Casimir of the jet (= jV c for gluons), and /i 2 /X g oc a 2 p is 
a transport coefficient of the medium proportional to the parton density, p. The 
factor, v ~ 1 — 3 depends logarithmically on L and the color Debye screening scale, 
/i. It is the radiated gluon mean free path, X g , that enters above. 

In the "thin plasma" approximation [ H [8] , the opacity expansion was applied 
and in the first order the following expression was derived for the energy loss: 



7T A, ./,, ki 

d V £ff 2k ± -q ± (k-q)lL 2 
Trfai+A- 2 ) 2 ' l6x 2 E 2 + (k-q)^^ 2 ' 1 ' 

Here the opacity factor L/X g is the average number of final state interactions that 
the radiated gluons suffer in the plasma. The upper transverse kinematic limit is 

k 2 nax = min [4E 2 x 2 , AE 2 x{l - x)] , (3) 

and the upper kinematic bound on the momentum transfer is <7n iax = s/4 ~ 3E[i. 
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Furthermore, A* 2 /// A* 2 = l+A'V'Zmax- For SPS and RHIC energies, these finite limits 
cannot be ignored [0 |5]. The integral averages over a screened Yukawa interaction 
with scale ji. The integrand is for an exponential density profile, p oc exp(— 2z/L), 
with the same mean thickness, L/2, as a uniform slab of plasma of width L. 

It was shown in Refs. [ [3 [8] that in the asymptotic E — > oo limit, the first- 
order expression ([2|) reduces to the BDMS result [ [3j S] up to a logarithmic factor 
log(E/fi). Numerical solutions revealed that second and third order in opacity 
corrections to eq. (j5J| remain small (< %20) in the kinematic range of interest. 

For applications to RHIC energies (y/s — 130 — 200 AGeV), we can find opacity 
values of L/X — 3 — 4 [ 113], which are connected to gluon rapidity densities of 
dN g /dy = 1000-1200 [[T4]. These values are consistent with the HIJING estimates 
on gluon rapidity densities [ [HI [16] and the measured charged hadron rapidity 
densities (~ 550 — 650) [[17l [18]. We illustrate here the jet quenching effects for 
a generic plasma with an average screening scale [i = 0.5 GeV, a s — 0.3, and an 
average gluon mean free path X g = 1 fm. The numerical results for the first order 
energy loss AE, taking into account the finite kinematic bounds are displayed in 
Fig. 1 for different opacities n = L/X g = 1 — 4. Since AE oc Cr, quark jet energy 
loss is simply 4/9 of the gluon energy loss shown in Fig. 1. As noted in Refs. [ 
[7] [5], the first order opacity result reproduces the characteristic BDMS quadratic 
dependence of the energy loss on L. 




E (GeV) 

Fig. 1. Absolute (AE) and relative (AE/E) energy loss of a gluon jet at different 
opacities, n — L/\ — 1, 2, 3, 4, as the function of the jet energy. 
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The most surprising feature of Fig. 1 is that unlike in the asymptotic BDMS 
case, where AE/E oc 1/E decreases rapidly with energy, the finite kinematic bounds 
in the GLV case give rise to an approximate linear energy dependence of AE in the 
energy range, E = 2 — 10 GeV. For this plasma, characterized by the transport 
coefficient fJ- 2 /X g = 0.25 GeV 2 /fm, the approximately constant fractional energy 
loss of gluon jets is 



^v/i-ljj^lljs) ■ (4) 



For the energy range displayed, AE B dms/E « (40GeV/i?)(L/6 fm) 2 for v = 2, 
exceeds unity throughout this energy range. The GLV expression approaches the 
asymptotic BDMS result from below only beyond the range of RHIC experiments. 

In order to investigate the influence of the GLV energy-dependent radiative 
energy loss on hadron production, we apply a perturbative QCD (pQCD) based 
description of Au + Au collisions, including energy loss prior to hadronization. First, 
we check that the applied pQCD description reproduces data on pion and kaon 
production in p + p collision. Our results are based on a leading order (LO) pQCD 
analysis. Detailed discussion of the formalism is published in Refs. [[Til H2]- Next 
to leading order calculations were performed for pions [ 1 13] . but not for kaons. 

Our pQCD calculations incorporate the parton transverse momentum ( "intrin- 
sic kr") via a Gaussian transverse momentum distribution g(kx) (characterized by 
the width (k T )) [HQ UJ as: 



b) x 



E ^~p~ J dx 1 dx 2 d 2 kT,ad 2 k T ,b g(k T , a )g(k T , 

abed 

ft r>2w i rfi\ da D h/c(z c ,Q 2 ) 

Ta/p{xi,Q )h/p{x 2 ,Q ) —j . (5) 

Here we use LO parton distribution functions (PDF) from the MRST98 param- 
eterization [[21] and a LO set of fragmentation functions (FF) [[22]. The applied 
scales are Q = p c /2 and Q = pt/2z c . 

Utilizing available high transverse-momentum (2 < px < 10 GeV/c) p + p data 
on pion and kaon production at 19 < ^fs < 63 GeV we can determine the best fitting 
energy-dependent (k 2 ,) parameter (see Ref. [[TTJ for further details on (fcf-)). Fig. 2 
shows a general agreement between the data and our calculations for the tt~ /n + 
and the K~ / K + ratios within the present errors for pt > 3 GeV. We also show 
predictions for y/s = 130 (56) GeV, with (k 2 ,) = 3.5 (3.0) GeV 2 /c 2 for tt+, tt" and 
K ~ . For K + we used a smaller phenomenological value, (k^) — 2.5 (2.0) GeV 2 /c 2 , 
as at lower energies [112]. 

Fig. 3 displays the results for K~ /ir~ and K + /tt + . One can conclude that, 
while agreement is reasonable for px > 4 GeV, there is a systematic discrepancy at 
smaller transverse momenta. High statistics p + p data at RHIC will be essential to 
establish an accurate baseline to which A + A must be compared. 
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Now we turn to the calculation for Au + Au collision at RHIC energy, y/s = 130 
AGeV. As a first approximation, let us consider slab geometry, neglecting radial 
dependence. We include the isospin asymmetry and the nuclear modification (shad- 
owing) into the nuclear PDF. We consider the average nuclear dependence of the 
PDF, and apply a scale independent parameterization with the shadowing function 
So,/a(x) taken from Ref. [US]: 



fa/A(x) = S a / A (x) 



z ( z 

■jfafp(x) + f 1 - -J ) fa/n(x) 



(6) 



The value of the (fcf.) of the transverse component of the PDF will be increased 
by multiscattering effects in A + A collisions. Two limiting cases were investigated 
with (i) a large number of rescatterings [ [12l [19] [20] and (ii) a small number of 
rescatterings ( "saturated Cronin effect" ) [ [11] . To clarify the situation, further 
study of p + A collisions is necessary. Here we concentrate on the influence of jet- 
quenching on the hadron spectra, leaving the inclusion of the multiscattering effect 
for future work. 

In this simplified calculation, jet quenching reduces the energy of the jet before 
fragmentation. We concentrate on y cm = 0, where the jet transverse momentum 
before fragmentation is shifted by the energy loss, p*(L/X) — p c — AE(E, L). This 
shifts the z c parameter in the integrand to z* = z c /(l — AE/p c ). The applied scale 
in the FF is similarly modified, Q = pt/2z*, while for the elementary hard reaction 
the scale remains Q = p c /2. 

With these approximations the invariant cross section of hadron production in 
central A + A collision is given by 

Eh J; = B { aa^Z Jdxidx2d 2 kT, a d 2 k T ,b g(kT,a)g(kT,b) 

abed 

t 1 ^>2w / n 2-s d(T '4 D h/c(z*,Q 2 ) 

fa/A{Xl,Q )f b /A{X2,Q —J , (7) 

dt Z c TTZ C 

where B^ A — 2tt J fcmox b db T AA (b) with b max = 3 fm for central Au + Au collision 
and T AA (b) is the thickness function. The factor z*/z c appears because of the in- 
medium modification of the fragmentation function [|23]. Thus, the invariant cross 
section ([7} will depend on the average opacity or collision number, n — L/X g . The 
calculated spectra for pions and kaons are displayed for n — 0, 2, 3, 4 in Fig. 4 (top 
panels), with their ratios to the non-quenched spectra (bottom panels). 

We note that the GLV energy-dependent energy loss leads to a rather struc- 
tureless downward shift of the single inclusive yields of all hadrons because AE/E 
is approximately constant. This is in distinction to the estimates in [[15], where an 
energy dependent fractional energy loss, AE/E = (0.5 GeV/-E)(L/fm), was used. 
Our new result indicates that over the entire accessible energy range, the GLV ra- 
diative energy loss reduces the pt distribution by up to an order of magnitude for 
pions and somewhat less for kaons. 
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Fig. 4. Top: 7r + and K + spectra in Au + Au collision at y/s = 130 AGeV without jet- 
quenching (dashed line) and with n = 2, 3, 4 (full lines). Bottom: Ratios of the spectra 
with quenching to the "non- quenching" case for n — 1, 2, 3, 4. For tt~ and K~ one obtains 
similar spectra and ratios. 

Fig. 4 indicates the importance of the absolutely normalized hadronic spectra 
to study the effect of jet-quenching. The effect on particle ratios is, however, more 
interesting, as shown by Fig.-s 5 and 6. 

While the influence on ir~ /ir + is very small, the K~ /K + ratio is predicted to 
drop dramatically in the few GeV domain. The shift of solid K~ / K + curves to 
smaller px (relative to the dotted curve) can be used to test the nonlinearity of the 
energy loss as a function of L. It also provides a constraint on the magnitude of the 
/Lt 2 /A ff transport coefficient. The most conspicuous quenching effect is limited in the 
K~ /ir~ ratio to below pr < 4— 5 GeV/c. There is about a factor of 2 enhancement 
of the K + /ir + ratio for opacity n — 4 throughout this range. 
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These patterns arise from the interplay of the energy-shifted gluon and quark 
fragmentation functions into pions and kaons. Since the quark to gluon jet ratio at 
fixed pr depends sensitively on the beam energy for s/s — 20 — 200 AGeV, these 
patterns are expected to shift in a systematic way with bombarding energy as well. 

In summary, we applied the GLV parton energy loss [ [3 [8] in LO pQCD to 
pion and kaon production in nuclear collisions and calculated the sensitivity of the 
hadronic ratios to this property of the QCD plasma. The K + /ir + and K~ /K + 
ratios in the 2 < pt < 10 GeV/c range were found to be the most sensitive to the 
gluon opacity at RHIC energies. 

We want to compare our theoretical results to experimental data at RHIC 
energies. Unfortunately, final data are available only at pr < 2 GeV/c in Au + Au 
collisions at yfs = 200 AGeV [ 24J . Looking into these data we find surprisingly large 
values for the kaon to pion ratios at px — 2 GeV/c in the most central collisions 
(0-5 %), namely K+/tt+ = 0.75 ± 0.05 and K~ /tt~ = 0.65 ± 0.05. Even in 
peripheral collisions (60 — 92 %) a moderate value, 0.45 ± 0.07 has been measured in 
both cases. These values are larger then any of our calculated ratios up to opacity 
L/X = 4 in the low momentum region, see Figure 6. This indicates that pion and 
kaon production is not dominated by fragmentation in this transverse momentum 
region. Furthermore, at lowpr we see from the data that K~ /K + « 1 [H3]. Figure 
5 displays that this ratio should drop to a small value around pr = 4 GeV/c in case 
of jet fragmentation. Thus one can determine from the drop of K~ /K + ratio the 
recovery of the fragmentation region. It is interesting to note that p/p ratio is 
constant (w 0.75) up to pr < 6 GeV/c [[24j [25J [26], which may indicate the lower 
limit of the fragmentation region. The measured anomalously high p/n + and p/ir~ 
ratios in the intermediate momentum region 2 < px < 6 GeV/c at RHIC energies [ 
125] have been explained theoretically by quark coalescence/recombination models [ 
[27] [28l [29] . Thus we may expect the application of similar models for mesonic 
ratios in the same momentum region. However, for pt > 6 GeV/c our calculations 
with fragmentation will be valid and applicable to extract the opacity of the hot 
dense parton matter produced in Au+Au collisions. 
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